Linear dark energy equation of state revealed by supernovae? 



Vincenzo Salzano 1 , Yun Wang 2 , Irene Sendra 1 , and Ruth Lazkoz 1 
1 Fisika Teorikoaren eta Zientziaren Historia Saila, Zientzia eta Teknologia Fakultatea, 
Euskal Herriko Unibertsitatea, 644 Posta Kutxatila, 48080 Bilbao, Spain and 
2 Homer L. Dodge Department of Physics & Astronomy, 
Univ. of Oklahoma, 440 W Brooks St., Norman, OK 73019, U.S.A. 
(Dated: November 7, 2012) 

In this letter we propose a test to detect the linearity of the dark energy equation of state, and 
apply it to two different Type la Supernova (SN la) data sets, Union2.1 and SNLS3. We find that: 
a. current SN la data are well described by a dark energy equation of state linear in the cosmic 
scale factor a, at least up to a redshift z = 1, independent of the pivot points chosen for the linear 
relation; b. there is no significant evidence of any deviation from linearity. This apparent linearity 
may reflect the limit of dark energy information extractable from current SN la data. 
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Introduction - After the discovery of the accelerated 
expansion of our Universe pQ, many models have been 
proposed to solve the mystery of dark energy, both in 
the context of general relativity and of alternative grav- 
ity theories [2j . Models are characterized by an analytical 
expression for the dark energy equation of state (EoS), 
derived either by phenomenological or theoretical consid- 
erations; and generally depending on two or more param- 
eters [3] which can be constrained using observations. A 
very short, not exhaustive, summary of the most used 
EoS parameterizations is in [4, 5j; among them, the CPL 
model 4J is considered as the reference model. It is de- 
fined as w(a) = Wq + (1 — a)w a : Wq is the EoS present 
value (a = 1 or z = 0) and Woo = wq + w a is the asymp- 
totic value of EoS at early times (a — > or z — > oo). This 
model has some well known problems: the high correla- 
tion between its two parameters, and the high-redshift 
dependence of the parameter w a , which makes its use 
with low-redshift limited data questionable. An alterna- 
tive parametrization is Wang's model [5J: 
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w{a) ={—l) W0+ {—l) Wc > (1) 

where w c = w(a c ) is the EoS calculated at the pivot 
value for the cosmic scale factor a, chosen to minimize the 
correlation between the EoS parameters, wq and w c . In 
a broader view, we can see that Wang's model is a linear 
interpolation between two points, {a\,wi = w(ai)) and 
(a2,u>2 = iu(d2)); thus we can define a General Linear 
EoS model (GL) by 
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We obtain Eq. ([!]) if a± = ao = 1 and a 2 = a c . 

GL models can be used to test whether the dark energy 
EoS is linear or not in the scale factor. This test has 
nontrivial implications: if linearity is verified, all models 
in [5] deviating from it may be discarded. If the EoS is 



not effectively linear, then the CPL model is intrinsically 
wrong. To test this possibility in the most general way, we 
can generalize the GL model to higher-order polynomial 
functions; an interpolation function between M points 
can be written in a general fashion as: 

M ( M \ 

with M = 2 for the GL model; N = 3 for a second- 
order interpolation model (2NL); M = 4 for a third-order 
interpolation model (3NL); and so on. A third possibility 
is that current SNela data do not allow us to distinguish 
between a linear and a nonlinear EoS, and linearity could 
be an artifact coming from the data used. 

This test can be implemented very easily, as illustrated 
in Fig.[l] Let us suppose to have an exact EoS whose (un- 
known) profile is the solid line in the figure; the dots are 
the pivot values where the GL EoS is to be calculated; 
the dashed lines represent the GL models (linear interpo- 
lations) using different pairs of pivots. If the underlying 
EoS is not linear (left panel) , then we will obtain different 
estimations of the pivot values when using different GL 
parameterizations; on the other side, if the EoS is linear 
(or with a very small curvature; right panel), such esti- 
mations will yield the same linear relation, independent 
of the GL model considered. 

Data - We will use SNela data from two different sam- 
ples: the Union2.1 compilation [7] and the SNLS3 com- 
pilation that includes the three year data from the Su- 
perNova Legacy Survey [S]. The x 2 is generally defined 
as x 2 = A:F • C- 1 • Kt. AF = F theo - T ohs is the 
difference between the observed and theoretical value of 
the observable quantity, T\ this will be the distance mod- 
ulus fi for Union2.1 and the SNela magnitude m mo d for 
SNLS3: 

A* = 51og 10 [i/ di(^, ^m! #)] + Mo ; (4) 
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Figure 1. Test of linearity of dark energy EoS through the GL models denned in Eq. |2j: solid line is the underlying unknown 
EoS function which we want to reconstruct; dots are the pivot scale factor values chosen to define various GL models; dashed 
lines are the GL models calculated varying the pivot values. In the left panel, the underlying EoS is not linear; in the right 
one, it is almost linear (very small curvature). 
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with the dark energy EoS parameters vector (0 = 
(u>i,u> 2 ) for GL models). 

The matter content will be fixed as Q m = 0.25!^] [5]; 
and M. are nuisance parameters combining the Hubble 
constant Hq and the absolute magnitude of a fiducial SN 
la, and we minimized x 2 (through a Monte Carlo Markov 
Chain (MCMC) algorithm) by marginalizing over them 
[5]. Finally, C is the covariance matrix: in the SNLS3 
case, it depends on the parameters a and /3, considered 
as free fitting parameters; in the Union2.1 case they are 
fixed. For a more detailed discussion about the difference 
between the two approaches, we refer to [8]. To test each 
model, we calculate the Bayesian evidence, defined as the 
probability of the data D given the model M with a set of 
parameters 0, £{M) = J d0L(D\0,M)ir(0\M): ir(0\M) 



1 We fix the value of fi m in order to reduce the degeneracy be- 
tween cosmological parameters and EoS dark energy parameters 
(wq, w a ). We have explored also the case of f! m = 0.271 and Q m 
free. We omit the description of such results: we have a small 
shift in the absolute values of the cosmological parameters in 
question and, when Q m is free, somewhat larger errors on them, 
but our conclusions are not changed at all. 



is the prior on the set of parameters, normalized to unity, 
and L(D\0, M) is the likelihood function. The evidence 
is estimated using the algorithm in [10 . To reduce the 
statistical noise we run the algorithm many times ob- 
taining a distribution of ~ 1000 values from which we 
extract the best value of the evidence as the mean of 
such distribution. Then, we calculate the Bayes Factor, 
defined as the ratio of evidences of two models, Mi and 
Mj, Bij = £i/£j. If Bij > 1, model Mi is preferred 
over Mj, given the data. The Bayesian evidence is in- 
terpreted using Jeffreys' Scale: if In Bij < 1 there is not 
significant preference for the model with the highest ev- 
idence; if 1 < In B^ < 2.5 the preference is substantial; 
if 2.5 < \nBij < 5 it is strong; if \nBij > 5 it is decisive. 
We will use the (0 — 0.5) GL model as reference model i. 

Results - First we focus on the information that can 
be extracted using different GL models. Final results are 
summarized in Table |T] We have chosen 5 pivot values 
for the scale factor, equally spaced in the redshift space 
(an arbitrary choice which has no effect on final results) 
and up to a maximum redshift z = 1 (the redshift value 
where the SNela number density is statistically signifi- 
cant). They are z = {0,0.25,0.5,0.75, 1} (corresponding 
to scale factors a = {1,0.8,0.67,0.57,0.5}). For each GL 
parametrization we have two primary parameters (bold 
text in tables), namely, the EoS calculated at the pivot 
redshifts and directly derived as fitting parameters; and 
secondary parameters (plain text), the EoS calculated in 
the remaining chosen pivot values and derived from the 
obtained best GL EoS relation. We also perform fits with 
a cosmological constant and the CPL model for compar- 
ison. 

We find the following about the absolute value of the 
EoS parameters: 
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Figure 2. Summary of all GL models: SNLS3 sample joining 
al the GL models in Table. |U 

• the primary EoS parameters are almost perfectly 
consistent with each other, independent of the GL 
model considered. We note that, given the nature 
of the MCMC algorithm, there is an intrinsic sta- 
tistical background noise (i.e. nonphysically mean- 
ingful fluctuations) to be considered; 

• the same good agreement is valid for secondary pa- 
rameters; particularly impressive is the consistency 
with the Woo = wq + w a derived from the CPL 
model. 

These are two strong points in favor of a linear (in scale 
factor) parametrization for the dark energy EoS, at least 
up to z = 1. Of course, we have to compare such results 
with the errors on these parameters: 

• we join all GL models: for each GL parametriza- 
tion, we randomly extract N = 1000 pairs of pivot 
EoS parameters, (u>i, 1^2)1 from Gaussian distri- 
butions centered on the best fit values from each 
model and with dispersion equal to the MCMC de- 
rived errors; then we plot the resulting GL models 
in Fig. [2] and we show the total errors on each EoS 
parameters in the last row of Table [I] We see there 
is a minimum in dark energy EoS parameters error 
at z ~ 0.25; this same value is obtained in [5] as 
the redshift corresponding to the minimal correla- 
tion between the pivot parameters, while errors on 
wo are comparable with errors on ^0.75; 

• the errors generally grow when moving to higher 
redshift pivot values: this is somewhat expected 
as, in general the SN la number density fades out 
toward higher redshifts. If we consider 0.1-width 
redshift bins, the SNLS3 sample has ~ 120 objects 
at z < 0.1 and ~ 40 objects per bin up to z ~ 
0.9; this, of course, influences the estimation of the 



pivot EoS parameters with a degradation of the 
reconstructed value at higher redshifts; 

Finally, we note that the Bayes factor (with the (0 — 0.5) 
GL model as reference) is much less than 0.05: GL mod- 
els are practically equivalent. This is an important con- 
sistency check in favor of linearity. But, considering the 
behavior of the errors, an important question arises: does 
the error- redshift relation depend on the number of data 
points, or is it a signal of an intrinsic breakdown of the 
linear EoS model?. As we have shown in the left panel of 
Fig. [TJ a non-linear EoS can produce different estimated 
values for the pivot parameters or, equivalently, a wider 
dispersion, i.e. larger errors on them. In order to check 
this, we have analyzed the 2NL and 3NL models, so that 
we can compare linear models to second and third-order 
interpolation functions. Results are given in Table. [TTJ 

• the primary EoS parameters are not consistent with 
each other as in the case of GL models; in particular 
the agreement among the derived values goes 
down; 

• errors are much larger than the GL models, due to 
the larger number of estimated parameters; 

• interestingly, even with these models, the W0.25 
pivot corresponds to the best estimation of the dark 
energy EoS; this could be a further indication that 
this depends on the higher number of points at 
lower redshift more than on intrinsic properties of 
EoS; 

• it is quite puzzling to note the value of the Bayes 
Factor: it is > 0.4 for the 2NL models and > 0.35 
for the 3NL ones. It is quite larger than the val- 
ues from the GL models, so that non-linear models 
are disfavored; but it is also lower than the val- 
ues which, following the Jeffreys scale, would give 
strong preference to the linear model. 

All the points above lead to the conclusion that a lin- 
ear EoS for dark energy is the most statistically proba- 
ble choice when using SNela as cosmological tools. The 
black region in Fig. [2] is the 68% confidence region, the 
most probable region where w(a) can lie. Clearly the val- 
ues of the errors, as explained above, still leave open the 
possibility for non-linear models; probably, with the ad- 
dition of more SNela in the future and at higher redshifts, 
we will be able to give even stronger constraints to the 
global EoS trend (presumably reducing the black region 
in the figure) and we would probably enhance further the 
evidence for linearity of dark energy EoS. 

Nevertheless, we have to consider another possibility. 
We have also used the Union2.1 sample; results for it 
are given in Table. |III| All the comments we have given 
for SNLS3 are qualitatively and globally true also for the 
Union2.1 data. In general, Union2.1 data show errors 
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which are always slightly smaller than those from SNLS3, 
mainly because Union2.1 has more low redshift SNela (if 
we consider 0.1-width redshift bins, we have ~ 100 ob- 
jects at z < 0.1, then we have a peak of ~ 60 objects 
per bin in 0.1 < z < 0.5 and ~ 20 objects per bin in 
0.7 < z < 1.0) which can more tightly constrain the 
behavior of EoS at higher redshifts. The behavior of 
the GL models at high redshift is even more interesting; 
we obtain two completely different asymptotical behav- 
iors: SNLS3 points to a very phantom early dark energy 
(woo ~ —2.6) which implies w a < and \w a \ > \wq\; 
Union2.1 has less negative EoS (w^ = —0.7), or equiv- 
alently, w a > and \w a \ < \wq\. While flux-averaging 
[TT] may reduce this dramatic difference somewhat, it 
will likely not remove it, since it most probably origi- 
nates from the fact that covariance matrices are obtained 
and handled very differently when obtaining the two data 
samples. However, both samples are very compatible 
with a linear EoS. 

To conclude, we have found that current SN la data 
are well described by a dark energy EoS linear in a, inde- 
pendent of the pivot points chosen for the linear relation, 
and that there is no significant evidence of any deviation 
from linearity. This may indicate that the dark energy 
EoS is a linear function in a, or that current data only 
allow the extraction of a linear function. Significantly 
larger SN la data sets will be required to clarify this. 
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Table I. Linear EoS test with SNLS3 data. For any EoS parametrization (first column) we give primary fitted parameters in 
bold text and secondary derived parameters in plain text. Errors for the primary parameters are derived directly from the 
MCMC procedure; the errors on secondary parameters are calculated following the standard error propagation theory. 
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Table II. Non-linear EoS test with SNLS3 data: column description is the same of Table. [I] 
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*(0 - 0.75 - 1.0) 


-0.983lS;|ii 


-1.150±1;11? 


-1.549iJ;g| 


-1 Q74+ 1 129 

X.IJ 1 1_2.142 


-2.370tJ;| 4 J 


-6.989tg;™ 7 


1 -tqo+0.109 

IA6S_ 105 


3.2861°;"? 


0.39 


(0.25 - 0.5 - 0.75) 


1 aoq+4.544 

±.uoo_ 6 g3S 


-1.168t°;=g 


— 1.560^ 788 


-1.9881^? 


-2.390t 2 3 ; 4 7 6 


-7.1461!°;™! 


± -^°^-0.100 


3.287; ;;;? 


0.49 


*(0.25 - 0.5 - 1.0) 


-0.852^j« 




loJ5 -0 .876 


-1.856±°£ 


-2.130tJ; 7 ^ 


-4 qn3+ 16 - 540 


1 4 or+0.113 
ll4Od -0 .104 


O OOQ + " 116 
°- zoo -0.111 


0.41 


(0.25 - 0.75 - 1.0) 


-1.0631S,, 


_1 -,7^ + 0.245 




-2.003li^ 


-2 41 ^ + 1 969 

Z ' 4i!5 -2 .947 


-7 Q49+ 41340 
' -°^— 58.128 


±.tO^_Q 102 


oQ/i+0.114 
o.z,o^_ .io7 


0.47 


(0.5 - 0.75 - 1.0) 


-i mn+ 31 - 507 

1.1UU_ 43 681 


-1 204+ 7 - 430 

l.ZU^_ 10 429 


-1 rne;+° 393 
±.ouo_ 62g 


-1.834il;°™ 


-2.143t 2 ;«K 


c onn+ 68 - 838 


l-434i' ;i^ 


3.280l° ;» 


0.38 


all (2NL) 


-1 014 +2 - 155 


1 isn+ - 635 

-l-loU_0 .627 


-1 604+ ' 865 


-2.036tl'j?= 


-2.45lt 2 ; 7 8 7 


-7 4QO+ 33 - 305 

' -^ yz -34.154 


- 


- 


0.46 


(0 - 0.25 - 0.5 - 0.75) 


-1 604+° 875 


1 9QQ+ - 300 


-2.0911™ 8 , 4 


1 2qq +1 193 




87 4qn+ lt)115e 

°'-^ yu -156.963 




oqo+0.113 
°- ouo -0.113 


0.35 


(0 - 0.25 - 0.5 - 1.0) 


1 cfto + 0.906 
i-0«0_o 99 6 


1 OTT+ 293 
L.OX i _ .34o 


-2.045;°,;^ 


-1.310il|S 


n ^7fi+ 3 841 


on 007+91. 018 
BU-OO' -109.622 


1 436+ 112 
± -^ ou -0.107 


3.303l ;" 4 


0.37 


*(0 - 0.25 - 0.75 - 1.0) 


-1384t°;l« 


— i.zui_ 316 


-1-8211KB 


• fi - LOO -6.167 


-2 1fi4 +4 730 

^ AD ^-12.359 


1 3 904+223.327 
iJ - zz4 -634.365 


1 407+0.112 


3.284; 4 i 9 


0.56 


(0 - 0.5 - 0.75 - 1.0) 


1 ^SQ+° 856 
1 ' UOB -0.903 


-1 4fifi+ H - 623 


^•UO1_ 900 


' ' °-1.864 


7fi7+ 3 055 
u - 1 ' -4.616 


77 8,1 +229.220 
' ' ■ ool -367.508 


1 494+OH2 

i-^-J^-o-ios 


3.303l ;5; 2 


0.25 


(0.25 - 0.5 - 0.75 - 1.0) 


1? 21 o+38.421 


1 qne+0294 

— i.^^o_ 301 




1 12 n+1114 
— 2.076 


-i.o78;| | T 


-100.955lg|; 2 ^ 


l.^oj_ 105 


oAo+0.117 
°- ouo -0.110 


0.34 


all (3NL) 


-1 436+ 1559 

i '^ JU -l .167 


-1 sm+ - 492 


-2 ofn +H " 7 


-1.388±igS 


-0.252l 4 ;f 23 


fi2 10S+ 194079 






0.37 



7 



Table III. Linear EoS test with Union2.1 data: column description is the same of Table. [I] 



Union2.1 


w 




W0.25 


Wo. 5 


Wo. 75 


Wl 


Woo 


CPL 


-0.933^ 


0.231 
0.224 


—0 qn7+ ' 327 

U.9U<_ 344 


-0.890t°; 4 I 5 2 


-0.877: 


r-0.592 
-0.642 


-0.8681°;?"! 


— 0.8031* |e 9 


(0-0.25) 


-0.972^ 


0.246 
0.246 


-0 912+° 105 

U - 3 - L ^ 0.113 


-0.8721S 


-0.843: 


H0.360 
-0.371 


-0.822±g;«f 


—0 fi72 +1 - 115 

u -°' z -1.135 


(0-0.5) 


-0.967^ 


0.245 
0.247 


-0 91 5+° 167 


-O.SSOl , 22 * 


-0.855: 


hO.298 
-0.308 


-0.8361°;^ 


-0 706+ 834 

u - ' UD -0.856 


(0 - 0.75) 


-0.968^ 


0.244 
0.252 


-0 Qio+ - 204 


u -°'°-0.286 


-0.848: 


f0.337 
-0.360 


-0.8281°;^ 


-0.6881°;^ 


(0-1.0) 


-0.968^ 


0.244 
0.249 


u.aio_ 233 


-0.878lg5g 


-0.850: 


hO.369 
-0.384 


-0.830l° ;* 28 6 


-0.692l°; 8 9 f 6 


(0.25 - 0.5) 


-0.988^ 


-0.413 
-0.448 


-0 911+° 102 
u - a J -- L 0.110 


-0.860l° 2 ^ 


-0.824: 


r-0.379 
-0.412 


-0 796+ ' 505 

U. 1 »O_ .548 


-0.605ll;I?l 


(0.25 - 0.75) 


-0.972^ 


-0.357 
-0.378 


-0.9131^ 


-0.874+S 


-0.846: 


f0.340 
-0.360 


-0.825l° ;f 7 I 


-0.67911;^ 


(0.25 - 1.0) 


-0.966^ 


-0.329 
-0.341 


-0 911+° 100 

0.111 


-o.874iS:5S 


-0.848: 


HO. 325 
-0.329 


-0.828™ 


-0.690±1;^| 


(0.5-0.75) 


-0.969tl-.gS 


U.91Z_q 7 g 3 


-0.874l° 222 


-0.847: 


0.335 
-0.357 


-0.827t°;«» 


— O.685I2J90 


(0.5 - 1.0) 


-0.961^ 


-1.071 
1.151 


-0.9091°;™ 


— n 87^+° 218 

U O ' °-0.241 


-0.850: 


r-0.260 
-0.276 


-0.8321°;"* 


-0.703±i;S« 


(0.75 - 1.0) 


-0.954^ 


-3.505 
-3.672 


-0.906+i;^ 


-0 S75+ - 976 

u -°'°-1.025 


-0.852: 


f0.338 
-0.357 


-O KlQec+0.431 
UOOO -0.448 


-O.7i6l4.3f9 


all 


-0.963j 


-0.483 
-0.483 


— n qi 9+°- 228 

U.91Z_ 238 


-0 877+ ' 268 

U.O ( / -0.261 


-0.840: 


r-0.336 
-0.348 


-0 832+ ' 467 
v.ooz,_ 451 


-0.687ll;f 2 



